On the valence-bond solid phase of the crossed-chain quantum spin model 
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Using a series expansion based on the flow-equation method we study the ground state energy 
and the elementary triplet excitations of a generalized model of crossed spin- 1/2 chains starting 
from the limit of decoupled quadrumers. The triplet dispersion is shown to be very sensitive to 
the inter-quadrumer frustration, exhibiting a line of almost complete localization as well as lines 
of quantum phase transitions limiting the stability of the valence-bond solid phase. In the vicinity 
of the checkerboard-point a finite window of exchange couplings is found with a non-zero spin-gap, 
consistent with known results from exact diagonalization. The ground state energy is lower than 
that of the bare quadrumer case for all exchange couplings investigated. In the limiting situation 
of the fully frustrated checkerboard magnet our results agree with earlier series expansion studies. 

PACS numbers: 75.10.Jm, 75.50.Ee, 75.40.-s 
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I. INTRODUCTION 

Frustrated quantum magnets have attracted consid- 
erable interest recently. On the one hand this is 
related to the ongoing quest for systems exhibiting 
spin liquid behavior and quantum disorder^. On the 
other hand numerous materials displaying geometrical 
or quantum chemical frustration of the magnetic ex- 
change have been discovered recently. Prominent ex- 
amples are the one-dimensional (ID) frustrated spin- 
Peierls compound CuGeOs^ or the 2D orthogonal spin- 
dimer system SrCu2 (603)2 with frustrating inter-dimer 
coupling^ as well as the 3D tetrahedral tellurate com- 
pounds Cu2Te205X2, with X=C1, Bri^. One promising 
route into spin-liquid behavior is via the coupling of lo- 
cally frustrated units like triangles or tetrahedra, as in 
the kagome, the pyrochlore, and the checkerboard (i.e. 
planar pyrochlore) lattices. Both, on the kagome and the 
checkerboard lattice low-lying singlets seem to exist with 
no long-range magnetic order (LRO) in the formeiiSii, and 
a valence-bond-crystal (VBC) ground state in the latter 
pg^gg8j9jio^_ The physics of the 3D pyrochlore quantum- 
magnet is still far from clarified with a hope that analysis 
of the planar model may lead to further progress in three 
dimcnsionsi2iiiiiii£. 

A closely related approach to spin-liquids above one 
dimensions has become of interest recently based on the 
frustrated coupling of spin-1/2 Heisenberg chains to form 
the planar crossed-chain model (CCM)iS*il'i^. A gener- 
alized version of this is depicted in fig. ^ which resem- 
bles the generic CCM for jo=ji- The generic CCM in- 
terpolates between three rather distinct regimes. Due 
to the complete frustration of the inter-chain exchange 

jo = ji E^nd for j2/ji ^ 1.25, it has been argued 
that the CCM stabilizes a 2D 'sliding Luttinger liq- 
uid' (SLL)iSiiLia. This is in contrast to the instabil- 
ity of coupled Heisenberg chains with respect to spinon- 
binding and the formation of antiferromagnetism (AFM) 
or VBC order if linked via norj-frustrating exchange. The 
SLL shows no LRO with massless, deconfined spinous 
forming the elementary excitations. At jo—ji^j2 the 
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FIG. 1: a) The generalized crossed-chain model. Spin-1/2 
moments are located on the solid circles. Thin solid(Thick 
solid and thin dashed) lines refer to the crossed chains(inter- 
chain coupling). The generic model with complete frustration 
of the inter chain coupling occurs at jo — ji . The quadrumer 
exchange is set to jo = 1 hereafter with ji and j2 correspond- 
ing to the expansion parameters, b) Local energy Ei, total 
spin Si, and quantum-number qi of a single quadrumer. 



CCM is identical to the checkerboard magnet exhibit- 
ing the aforementioned VBC ground state and a spin 
gapSiSiiSiii. Finally, for j2 and jo=ji the CCM maps 
onto the 2D spin-1/2 Heisenberg model on the square- 
lattice with AFM LRO and gapless magnon excitations. 
The AFM LRO in the latter region and the VBC have 
been suggested to coexist in the vicinity of the point 
j2/ji ^ 0.7^^'^", which however has been questioned 
based on results from exact diagonalization (ED)ii. 

The CCM is under intense current investigation. The 
SLL limit has been analyzed by bosonizationiiSii2ii&. 
Numerical studies have been performed on up to 36 
gpjj^g8j9jiij2^^ In particular, ED for jo=ji in ref>ii sug- 
gests the VBC spin-gap to close below j2/ji ^ 0.65 
and moreover ID, i.e. SLL, behavior is found above 
j2/ii ^ 1-5. Several analytical methods have been 
employed, including various semiclassical, Sp(N), and 
linear-spin- wave (LSW) approaches incorporating 1/S- 
corrections alsoiSiSSiSiSiiS^. LSW predicts stability of 
the AFM LRO for J2/J1 < 0.76 at jo=ji and spin-1/221. 
Because it is not even obvious that extrapolations from 
the large-S" or -N limit will lead to conclusive answers 
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for the quantum case, it seems highly desirable to ob- 
tain results from other analytic approaches applicable to 
5=1/2, such as eg. series expansions (SE). At ji—j2 SE 
has been performed starting from the quadrumer^ and 
the tetrahedral^^ limit. However, SE for ji ^ j2 is not 
available. 

In this context it is the purpose of this work to shed 
more light onto the quantum CCM by studying its el- 
ementary excitations, the stability of the VBC phase, 
and the ground state properties. This will be done by 
series expansions starting from the limit of decoupled 
quadrumers. The paper is organized as follows. In sec- 
tion two we describe our method of calculation. In sec- 
tion three and four we discuss the triplet excitations and 
the instabilities of the VBC phase with respect to triplet 
softening. Section five is devoted to the ground state en- 
ergy. A summary and some technical details are provided 
in the conclusions and the appendix. 



II. SERIES EXPANSION BY CONTINUOUS 
UNITARY TRANSFORMATION 

Recently it has been shown, that the checkerboard 
point of the CCM, i.e. ji~j2=^, can be described by SE 
in terms of the single coupling constant j—ji—j2, start- 
ing from the limit of decoupled quadrumes^. The SE was 
found to adiabatically connect the bare VBC at j = to 
a renormalized one at j = 1 without being hindered by 
quantum phase transitions. This is consistent with other 
findings of a VBC at ji=j2= l?'^?i" . Motivated by this 
we decompose the Hamiltonian of the CCM in a form 
adapted to a VBC symmetry breaking and normalized 
to the overall unit of energy jo which will be set equal to 
unity hereafter 



SiiSi+ii + ji(S2iSii+y + S31S41-1 



S31S2I + X + S41S1I + X) + j2(S2lS41 + y + SslSll + y 
S31S1I + X + S41S2I + X)] 

N 

^ Ho+ ^ (jlTi„ -t- j2T2n) 



(1) 

(2) 



Si\ refers to spin- 1/2 operators residing on the vertices 
I = 1 ... 4 of the quadrumer at site 1 with periodic bound- 
ary conditions on i, c.f. fig. ^). Hq is the sum over lo- 
cal quadrumer Hamiltonians the spectrum of which, c.f. 
fig. ^), consists of four equidistant levels which can be 
labeled by the local total spin and a number of lo- 
cal energy quanta qi. Hq displays an equidistant ladder 
spectrum labeled by Q = The Q=0 sector is the 

unperturbed ground state |0) of Hq, which is a VBC of 
quadrumer-singlets. The (5=l-sector of Hq consists of 
linear combinations of local S\=l single-particle excita- 
tions of the VBC with 91=1, where 1 runs over the lattice. 
At Q=2 the spectrum of Hq has total spin 5=0,1, or 2 



and is of multiparticle nature. Eg., for total 5=0 at Q=2 
it comprises of one-particle singlets with qi=2 and two- 
particle singlets constructed from triplets with qi=qm=l 
and 1 ^ m. In turn, the perturbing terms in J^l ji.2, 
can be understood as a sum of operators Ti„_2n which 
nonlocally create(destroy) n > (n < 0) quanta within 
the ladder spectrum of Hq. 

It has been shownSi242LSS42S that problems of type 
allow for perturbative analysis using a continuous unitary 
transformation generated by the flow equation method of 
Wegneia2fl. The unitarily rotated effective Hamiltonian 
Hcs readsSLSa 

H,s = Ho + Y^ J2 C{ni)W,n,W^,...Wm^ (3) 



n=l 



m| — n 
Af(m)=0 



where m = (mi . . .to„) with |m| = n is an n-tuple of 
integers, each in a range of m.i € {0, ±1, . . . , ±7V} and 
Wn — jiTin + j2T2n- In Contrast to of (^l, H^s con- 
serves the total number of quanta Q. This is evident from 
the constraint M(m) = ~ The amplitudes 

C (m) are rational numbers computed from the flow equa- 
tion method^^'^^. Explicit tabulationSi of the Tin shows 
that for the Hamiltonian in Q = 4. In the context of 
spin systems the flow equation approach has been applied 
successfully to ID and 2D dimer-models^SiS, as well as 
to 2D and 3D quadrumer-models^ 

Q-conservation of Hcs leads to a ground state en- 
ergy of Eg = (O|i/oft|0). Q-conservation also guarantees 
the Q=l-triplets to remain genuine one-particle states, 
i.e., their dispersion can be calculated via Ef^{k) = 



n,im = {ii,lm\Hes\n,00) - 
SimfioEf"^ are hopping matrix elements from site (0, 0) to 
site (Z, to) for a quadrumer excitation /i inserted into the 
unperturbed ground state, t^^im has to be evaluated on 
clusters with open boundary conditions large enough to 
embed all linked paths of length n connecting sites (0, 0) 
to (Z,to) at 0(n) of the perturbation. Ef'' = (0|i7off|0) 
on the t^^oo-cluster. States from sectors with Q > 1 will 
not only disperse via HcS but require the solution of an 
interacting problem. 



III. ELEMENTARY TRIPLET EXCITATIONS 

We have evaluated the one particle hopping matrix 
elements up to fifth order. They are listed explicitly 
in the appendix. Several of the resulting dispersions 
-BT(kjJi,i2) are shown in fig. [21 both, as a function of 
the wave vector along the irreducible wedge of the Bril- 
louin zone of the 2D square lattice and for various values 
of the coupling constants ji and j2- 

Two effects can be extracted from fig. |2 First, it 
demonstrates that the CCM 'localizes' triplet excitations 
on the line of largest inter-quadrumer frustration, i.e. for 
ji = j2- This can be observed both, from the three pan- 
els a), b), and c) where the dispersions seem completely 
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flat for ji=j2, but also from the explicit expressions for 
tijn in appendix IBl I.e., up to 0(4) and for I ^ m all 
hopping matrix elements are proportional to powers of 
(ji ~ 32)- For I ^ m and Im 7^ 10 this remains true even 
up to 0(5), where however tip ^ 0. Therefore the triplet- 
localization at maximum frustration is not complete and 
a weak hopping remains beyond 0(4). Inserting ji=j2 
into (|B1IB13| ) we find that -ET(k, ji = j^) is exactly iden- 
tical up to 0(5) with an 0(7) SE which is available for 
this case^. 

The second fact to note from fig. [21 is that the VBC 
ground-state is unstable with respect to triplet-softening, 
i.e. magnetic ordering at suitably chosen values of the 
coupHng constants. Depending on the ratio of 32/ ji > 
1 or < 1 this instability will occur either at a critical 
wave vector of —{tTjTt) or (0,0). Consistently with 
this the structure of exchange for the CCM, although 
topologically distinct, is bipartite both, for ji ^ 0, j2 ^ 
and for ji — > and J2 7^ 0. The instabilities arc studied 
in more detail in the next section. 

Finally, panel a) and b) of fig. [21 contain 0(4) SEs for a 
comparison with the 0(5) results. The difference is small 
only. For ji and j2 along path c) the difference is hardly 
visible and the 0(4) SE has not been displayed in the 
corresponding panel. The case of j/i = 0.5 and j2 ~ 
in panel c) is identical to a particular realization of the 
plaquette square-lattice Heisenberg magnet. The triplet 
dispersion of this system has been investigated with a 
rather different type of fifth-order plaquette SE by Koga 
and Kawakamj2-. Their result for the triplet dispersion, 
is included by the dots in panel c) and is in excellent 
agreement with our findings. 



IV. QUANTUM PHASE TRANSITIONS 

In the following we consider the lines of quantum phase 
transitions resulting from the closing of the elementary 
triplet-gap. This provides for an upper bound on the ex- 
tent of the region of stability of the VBC in a tentative 
quantum phase diagram of the CCM. In principle addi- 
tional first order transitions could occur or excited states 
other than the elementary triplets could collapse onto the 
ground state as well, leading to additional phase bound- 
aries. Here we focus on the triplet instability. Moreover, 
we limit our discussion to AFM couplings Ji(2) > and 
consider ferromagnetic cases elsewhere'^^. 

Before proceeding we stress that the CCM is symmetric 
under the interchange ji ^ jo and a diagonal shift by one 
unit cell. This implies that the complete quantum phase 
diagram can be obtained from a mirror reflection of the 
region < ji < 1 and < j2 at the line (ji — l,j2) and 
a subsequent relabeling of the axis by ji Jo/ji = 1/ji 
and j2 32/31- 

The results are summarized in figs. [3 and 01 which dis- 
play information from the plain series as well as from a 
Dlog-Pade analysis. In fig. [3^) we show the evolution 
with increasing order of the lines of vanishing triplet gap 
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FIG. 2: Evolution of the triplet dispersion i5T(k, 71,72) for 
various ji, 72 chosen along the directions a), b), and c) 
sketched in the top panel and for a path of k=(0, 0)-(7r, 0)- 
(tt, 7r)-(0, 0) in the 2D Brillouin zone. Solid[dashed] curves in 
panels a)-c) [a), b)]: 0(5)[0(4)] SEs. Dots in panel c): differ- 
ent 0(5) SE for the plaquette square-lattice from ref^. 



of the hare series. First, this figure further clarifies the 
location of the critical wave-vector kc as anticipated al- 
ready in fig. [21 Second the panel is intended to prove 
a monotonous and smooth behavior of the bare series 
within the range of interest. Based on the 0(5) SE, 
fig. shows the transition lines resulting from single- 
variable Dlog-Pades which have been obtained by pa- 
rameterizing j2(i) of the bare SE by 72(1) = aii(2) for 
j2(i) < ii(2)- The figure displays a [2,2] and a [1,3] Dlog- 
Pade. They are indistinguishable on the scale of the plot. 
The lower(upper) critical value of ^2 which they yield for 
the closing of the triplet gap on the line (ji=l,j2) are 
= 0.81 and = 1.06. As discussed in the next para- 
graph, along the diagonal, i.e. for ji=j2, Dlog-Pades can 
be obtained also from an 0(7) SE^. They result in some- 



4 




J 0.5 



J 0.6 



FIG. 3: Triplet instabilities of the generalized CCM: a) Lines 
of vanishing triplet-gap for increasing orders of the bare SE 
from 0(1) to 0(5). Triplet softening occurs at the wave vec- 
tor (tt, tt) and (0, 0) for the upper and lower set of lines. VBC 
labels the region of a finite triplet gap. b) Comparison of the 
fines of vanishing triplet-gap at 0(5) of the bare SE (solid) 
with the corresponding [2,2] and [1,3] Dlog-Pade approxi- 
mants (dashed-dotted). The latter two are indistinguishable 
on the scale of the plot. Along the selected directions 1), 
2) and 3) cuts from the Dlog-Pade analysis are be exempli- 
fied in fig. 13] The solid dot at ii=j2 = l refers the checker- 
board point. The regions labeled by italic style AFM{') are 
likely to exhibit antiferromagnetic ordering— On 
the line (ji — 1,72) and for j2 > the ground state is a 
'sliding Luttinger liquid (SLL). The remainder of this 
diagram for all coupling constants is obtained by reflection 
along [ji = 1,72) and a rescaling of ji — > and j2 32/31 
(see text). 



what larger values of the critical couplings than the 0(5) 
SE. Roughly interpolating the difference in between 
the 0(5) and 0(7) Dlog-Pades along this direction down 
to the line (ji=l,j2) we approximate an error for the criti- 
cal couplings of = 0.79 . . . 0.81 and = 1.06 . . . 1.13. 
These values should be contrasted against LSW theory^^, 
which predicts stability of an AFM phase for j2/ ji < 0.76 
at ji = 1. Moreover, ED at ji = 1 for up to 36 spins^^ 
suggests a closing of the spin-gap below j2/ji ~ 0.65 and 
a crossover to ID SLL-behavior at j2/ji ^ 1.5. 

The role of the Pade analysis in fixing the location of 
the quantum critical points as compared to the bare SE is 
clarified in more detail in fig.Qlwhere several reintegrated 
Dlog-Pades are displayed along with the plain series for 
three selected directions shown in fig. Eb). For ji = j2, 
where the Pade analysis is expected to be most relevant, 
an 0(7) SE from ref.^ has been included in fig.E]3). Be- 
cause of the symmetries of the quantum phase diagram 
the curves beyond the vertical dashed line in fig. E]3) do 
not refer to actual values of the spin gap for j > 1. This 
range is shown for completeness sake only. Most impor- 
tant, from fig. ^3) one realizes that both, the plain SE 
as well as all of the Pade approximants yield no critical 
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[3,1]-[1,3],0(5) 
[5,1]-[1,5],0(7) 



J 



FIG. 4: Spin gap At of the bare series (thick solid and 
thick dashed-dotted) vs. reintegrated Dlog-Pade approxi- 
mants along three selected directions 1), 2), and 3) shown 
in fig. 133), i.e. ji = j and j2 = aj with a = 0,0.5 and 1. 
The 0(7) SE in panel 3) is from refi^. At 0(5) all ( the [2,2] 
and [1,3] ) approximants are indistinguishable on the scale of 
panel 1) ( 2) and 3) ). At 0(5), in panel 2) and 3) (0(7), in 
panel 3) ) the [3,1] ([5,1]) approximant is split off from the 
remaining ones and is closer to the bare SE. 



points for j < 1. While this agrees with earlier find- 
ings of a finite spin gap ai j — 1 in refs.^'^'"'^^ it seems 
that the checkerboard point is almost critical. Panel 3) 
shows that the [4,2]-[l,5] and the [2,2]-[l,3] approximants 
at 0(7) and 0(5), respectively, are well clustered with 
critical points clearly separated from those of the corre- 
sponding plain SE and the [5,1] and [3,1] approximants. 
Therefore, the Dlog-Pade analysis at 0(5) and arbitrary 
ji, j2 in fig. Ob) iias been based on the [2,2] and [1,3] 
approximants. The effects of the Dlog-Pades in figs.^1) 
and 2) is rather minute, with the [2,2] and [1,3] approxi- 
mants leading to essentially identical spin gaps. 

As noted already, direction 1) in fig. Ob) refers to the 
plaquette square-lattice Heisenberg model. Its transition 
to the AFM state for j ^ 1 has been studied by several 
groups. We may therefore compare the critical coupling 
jc = 0.555 which we find with that obtained by other 
SE investigations, i.e. jc = 0.555 in ref.'^^ and jc — 0.54 
in refi^, as well as by ED, i.e. jc = 0.6 in ref.'^^. The 
agreement is satisfying. 

Regarding the symmetry of the ground state in region 
AFM(') in fig. Ob); rio conclusions can be drawn from 
the quadrumer SE. Yet, subsuming the present work with 
j.gfg Jj^22j35,36j37 jg ygj.y i[]^Qiy t}iat region AFM exhibits 

simple Neel-type AFM-LRO. Within AFM' the CCM ac- 
quires a bipartite lattice-structure as ji — *■ which could 
trigger a '4-up-4-down'-type AFM-LRO of the crossed 
chains for j2 > 1. The crossover between this and the 
SLL phase as ji — > 1 is an open issue. 
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FIG. 5: 0(6) ground-state energy eg per spin from eqn. JIJ. 
eg resembles the ground-state energy of the CCM only within 
the range of stability of the VBC bounded by thick solid lines 
in the (ji j2)-plane, i.e. the dashed-dotted lines from fig.Ot). 
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FIG. 6: Ground state energy of the CCM on the line (ji = 
1,^2). Thick solid, dashed-dotted, and dashed-double-dotted 
lines: SE. Solid curves with symbols: ED from fig. 2 of ref^. 
jj^'^* critical couplings from Dlog-Pade at 0(5), see fig.EJs). 



V. GROUND STATE ENERGY 

To obtain SEs for the ground state energy valid to Oin) 
in the thermodynamic limit, i.e. for systems of infinite 
size one may evaluate the matrix element (0|iJoff|0) on 
single 'maximal' clusters, i.e. with periodic boundary 
conditions and sufficiently large not to allow for wrap 
around at graph-length n^iSS. Here, Q-conservation of 
i/oft is the main ingredient which reduces the size of the 
intermediate Hilbert spaces. However, for the CCM we 
find the maximal-cluster approach to be severely limited 
by memory-constraints and we have combined the flow- 
equation approach with a linked cluster expansion. Some 
details of the latter are commented on in appendix^ Up 
to 0(6) the ground state energy per s'pin is 

1 481 Y19ih\j\ 5705977jt 

^ "2 ^ "576" " 13824 " 663552 ~ 522547200 
13033565594599jg 2nn ^^jjh 224267j3j2 
2123765592883200 9 1728 4354560 
388714973jfj2 655905584767ii^j2 67j| 

^ 35115171840 ^ 130026464870400 ~ ^76 
59jij| _ 4837781ji^j| 75848383^3^1 
1728 104509440 175575859200 
169095132323jtj1 481 j| 574241 

~ 168552824832000 ~ 13824 ^ 13063680 
38923349j?j| 113410023666229ji3j3 I39031j| 

58525286400 ^ 15928241946624000 ~ 5971968 
1184653457jij| 65325840449533ji2j4 2624063j| 

175575859200 ~ 10618827964416000 ~ 313528320 
181444615182347jij| 51497150603363j| 
^ 31856483893248000 ^ 10618827964416000 ^' 

First we note, that at the checkerboard point, i.e. for 
ji = J2 6qn. is identical to the corresponding SE for 
Cg in eqn. (6) of ref:-. Figure |31 visualizes e^, resem- 



bling a convex function with the weakest energy gain in 
the vicinity of the line of strongest frustration, i.e. for 
ji = 32- Note, that fig. [S] is the ground state energy of 
the CCM only within the region of stability of the VBC. 
Obviously eg is below the ground state energy of the bare 
quadrumer limit for all ji, j2. This is consistent with 
j-gfgj8j9jii^ and is at variance with a 'tetrahedral' SE to 
0(3) at ji=j2 in ref.^^, where Cg > —0.5 has been found. 

In rcf.^ it has been shown, that eg at ji—j2 agrees to 
within 1% with that of ED on systems with 36 spins for 
ji < 1 • In fig. El we compare eg with ED-results along 
the line (ji = 1,^2) obtained by Sindzingre, Fouet, and 
Lhuillier^ for systems with 16, 32, and 36 spins. The SE 
is displayed along with its variation in going from 0(4) 
to 0(6) in order to provide a rough estimate of conver- 
gence. As with fig. [SI the SE resembles the actual ground 
state energy only within the range of critical couplings 
^2^'^' marked on the j2-a.xis. While the line (ji=l,j2) is 
rather remote from the limit of decoupled quadrumers, 
the agreement between SE and ED on the largest system 
depicted in fig. El is still satisfying - even below j!^^ . 

Finally, one may define a 'point of maximum frustra- 
tion', i.e. j™^ at maximum Cg for ji = l. Interestingly j™^ 
^ 1, i.e. it is from the checkerboard point. We find 
that j^'^ « 0.973 ± 0.004. The error has been approxi- 
mated from twice the difference between j™''^ at 0(5) and 
0(6). As similar observation has been reported in ref^. 



VI. CONCLUSION 

To summarize, we have investigated the ground state 
properties and the elementary triplet excitations of 
the generalized crossed-chain quantum-spin model by 
quadrumer series expansion. Over a large region of 
exchange parameters the model exhibits a finite spin 



gap which is smoothly connected to the hmit of the 
bare quadrumers. This region is bounded by lines of 
triplet softening which results from a delocalization of 
the triplets away from the case of complete frustration 
of the inter-quadrumer exchange. This remains true in 
particular not only at - but also in the vicinity of the 
checkerboard point for ji=l where we have established 
a finite range of on-chain exchange, i.e. 0.79...0.81< 
j2 <1. 06. ..1.13 with a non-zero spin gap. While the fate 
of more complex multiparticle excitations as a function 
of ji(2) remains to be investigated in the future, our find- 
ings are consistent with a ground-state symmetry in the 
region of non-zero spin gap identical to that at ji(2)=0, 
i.e. a valence-bond-solid. 
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FIG. 7: Non-isomorphic (isomorphic) graphs of the 
quadrumer (standard) expansion on the crossed-chain (sim- 
ple square) lattice. Graph 2 is meant to results from 1 by: a) 
an in-plane 90° rotation of the upper edge, b) a 180°-out-off- 
plane rotation of the upper edge around edge 'e'. 



Only if ji=j2 cqns. lAll and IA2I are identical. In con- 
clusion, the number of non-isomorphic graphs in the 
quadrumer expansion for the CCM is significantly larger 
than that in standard expansions for 2D square-lattice 
models. 



APPENDIX A: LINKED CLUSTER EXPANSION 



APPENDIX B: TRIPLET HOPPING 
AMPLITUDES 



Here we comment on a peculiarity of graph counting for 
the quadrumer expansion for the CCM as compared to 
that for standard expansions for spin lattice-models with 
no internal structure of the vertices and edges. For the 
latter, both pairs of graphs in fig. [7^) and b) are isomor- 
phic and yield identical results. In general this is not the 
case for the quadrumer expansion. This is due to the ver- 
tices (edges) to consist of 4-spin quadrumers (tetrahedral 
links) which prohibit 'free rotation' of a graph at its ver- 
tices. Similarly, 'twisting' of an edge - as for edge 'e' in 
fig. {7}p) - is no symmetry operation in general. In fact, 
this depends on the choice of parameters, i.e. exactly at 
the checkerboard point, for ji=j2, the tetrahedral link is 
symmetric under the twist. This is refiected in, the cor- 
responding cluster weights for the ground state energies, 
which (discarding lattice constants) read 
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In this appendix, and for completeness sake, we list the 
triplet dispersion £^T(k) in the thermodynamic limit up 
to 0(j™, ^2 ) with m + n < 5 which can be written as 

EriKh^h) = ^00+ 

tio(cos{kx) + cos{ky)) + 

til cos{kx) cos{ky) + 

i2o(cos(2/ca;) -t- cos{2ky)) + 
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t^i cos{kx) cos(/cy)(cos(2fca;) + cos(2fcy) — 1) + 

t32{cos{3kx) cos{2ky) + cos(2A:a;) cos(3/cj,)) -|- 

t4o{cos{4:kx) + cos{4:ky)) + 

t4i{cos{4:kx) cos{ky) + cos{kx) cos(4fcy)) + 

tr,o{cos{5kx) + cos{5ky)) (Bl) 

where the hopping amplitudes ioo,...,50 are given by 
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